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Refractive Index Effects on 
Radiation in an Absorbing, 
Emitting, and Scattering 
Laminated Layer 

A simple set of equations is derived for predicting temperature distributions and 
radiative energy flow in a two-region semitransparent laminated layer in the limit 
of zero heat conduction . The composite is heated on its two sides by unequal amounts 
of incident radiation. The two layers of the composite have different refractive 
indices , and each material absorbs , emits , and isotropically scatters radiation. The 
interfaces are diffuse, and all interface reflections are included. To illustrate the 
thermal behavior that is readily calculated from the equations, typical results are 
given for various optical thicknesses and refractive indices of the layers. Internal 
reflections have a substantial effect on the temperature distribution and radiative 
heat flow. 


Introduction 

The use of ceramic materials for parts or coatings is of 
interest for high-temperature applications. Some of these ma- 
terials have reinforcing fibers or are laminated so it is necessary 
to consider heat transfer in composite regions. The surrounding 
temperatures are high, which provides substantial radiative 
heating. Since some of the materials are semitransparent, their 
temperature distributions depend on the internal radiative heat 
flow. The refractive index of a semitransparent material can 
have a considerable effect on its temperature distribution. The 
refractive index governs the amount of transmission into the 
interior of the material and affects the internal reflections that 
occur. Of more significance, the emission within a material 
depends on the square of its refractive index; hence internal 
emission can be many times that for a blackbody radiating 
into a vacuum. Since radiation exiting through an interface 
into a"vacuum cannot exceed that of a blackbody, the amount 
of internal reflection can be substantial. It redistributes energy 
within the layer and tends to make the temperature distribution 
more uniform than for a material with refractive index close 
to one. 

To obtain the general solution for the temperature distri- 
bution in a composite semitransparent layer requires solving 
the integral equations of energy transfer in each region, in- 
cluding heat conduction, and matching temperatures and heat 
flows at the internal interface. A spectral calculation can be 
carried out in each significant wavelength band and the total 
energy flow found by summing over the bands. A numerical 
solution by computer is required. A simple limiting solution 
is obtained here that is helpful and informative. The purpose 
of this paper is to show that in the limit with no heat conduction 
and for gray layers with refractive indices greater than one, 
the solution for a composite region can be obtained from 
information that is already known and results calculated very 
easily. The effects of surface reflections and isotropic scattering 
are included. The solution of integral equations is not required. 
The limiting result can be helpful in initiating iterative com- 
puter solutions for a more general analysis and for checking 
the validity of numerical solutions. 

The development here builds on the analysis of Siegel and 
Spuckler (1992) where it was shown that for radiative equilib- 
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rium in a gray layer with diffuse interfaces, the temperature 
distribution and radiative heat flux for any index of refraction 
is obtained very simply from available results for an index of 
refraction of one. The present paper applies these ideas to a 
two-layer laminate subjected to external heating on its outer 
surfaces. Each of the layers emits, absorbs, and isotropically 
scatters radiation. For simplicity the medium surrounding the 
laminate has a refractive index of one. 

The outer surfaces of the laminate, and the interface between 
the two layers, are each assumed diffuse. This is probably a 
reasonable approximation for unpolished materials that are 
bonded together. Transmitted radiation, or radiation emitted 
from the interior, is assumed to be diffuse when it reaches the 
inner surface of an interface. If the index of refraction of the 
material is greater than that of the surrounding medium, some 
of the internal radiaTion is in angular directions for which there 
is total internal reflection. This retains energy within the layer 
and tends to equalize its temperature distribution. 

Some of the important early work on calculating heat trans- 
fer behavior in semitransparent layers was done by Gardon 
(1958) who developed an analysis for temperature distributions 
during heat treating and cooling of glass plates including index 
of refraction effects. The interfaces were optically smooth, 
and the specular reflections at the interfaces were computed 
from the Fresnel reflection relations. A similar application 
(Fowle et al., 1969), predicted heating in a window of a re- 
entry vehicle. Recent papers by Rokhsaz and Dougherty (1989), 
Ping and Lallemand (1989), and Crosbie and Shieh (1990) 
further examined the effects of Fresnel boundary reflections 
and having a refractive index larger than one. Many analyses 
of both steady and transient heat transfer in single or multiple 
plane layers such as Amlin and Korpela (1979), and Tarshis 
et al. (1969), have used diffuse conditions at the interfaces as 
in the present study. Thomas (1992) set up a solution procedure 
to include a ceramic interface that is partially specular and 
diffuse. Ho and Ozisik (1987) carried out a numerical analysis 
of radiation and conduction in a two-region laminate without 
including refractive index effects and hence without internal 
interface reflections. 

Analysis 

A laminated plane layer is made of two different materials 
that have thicknesses D\ and D 2 as shown in Fig. 1. Both 
sublayers absorb, emit, and isotropically scatter radiation. The 
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Fig. 1 Laminated two-layer geometry, coordinate system, and nomen- 
clature for heat fluxes at interfaces 


limiting case is considered here where the temperature distri- 
bution within the composite is dominated by radiation so heat 
conduction is neglected. Each material has a constant refractive 
index; it is the effect of the two different n values that is 
investigated here along with having different optical thick- 
nesses for the layers. The materials provide significant scat- 
tering, so the interfaces between the two layers and between 
each layer and the surrounding air or vacuum (with n = 1) 
are assumed diffuse. The laminated layer is subjected to diffuse 
radiation from the surroundings of q° r \ and q° r2 on the two outer 
boundaries x x = 0 and x 2 = D 2 ; for convenience here q° n > 
q°r 2 . Inside each layer there are outgoing and incoming fluxes, 
q 0 and q h at each interior interface. Since scattering is included, 
the local optical depth in each layer is related to its x coordinate 
by k = (a + o s )x. Using an individual coordinate in each of 
the two layers facilitates the formulation and solution that 
follows. 

The temperature distribution in each layer is governed by 
an integral equation given by Siegel and Howell (1981), mod- 
ified with index of refraction factors as, 


rt?a 7 t(K|) = ^ [Qo.aEiM +q 0 ,bEi(K-D\-K\)] 

n 2 C* Dl 

+ y) o a7l ( k ,*)£,( Ik, - k‘ I )dK* (Iff) 

n\oT\ ( « 2 ) = 2 [<70 , <£i ( « 2 ) + Qo.dE 2 (k d1 - k,)] 
n 2 C* D2 

+ ~l\ °^ K2 )£l(l*2“K2 I )^2 GW 


The radiative flux is a constant through the entire laminated 
layer for the present conditions of radiative equilibrium. The 
relation of flux to temperature distribution in each of the two 
layers is obtained by evaluating the equation for q r for a plane 
layer from Siegel and Howell (1981) at x x = 0 and x 2 = 0, 

J K D1 

aT\{^K\)E 2 {K\)dK\ (2a) 

0 


J k D2 

oT 2 (K 2 )E 2 (K 2 )dK 2 (2b) 

0 


The following dimensionless groups are now defined: 


x n]oT\(K X ) - q 0tb ^ x n 2 2 aTt(K 2 ) - q oM LX 

*i(Ki) = ; ^ta) = — b) 

Qo,a ~~ Qo,b Qo,c Qo,d 


% 

Qo,a Qo,b Qo,c Qo,d 


(4 a, b) 


Equations (1) and (2) then transform to (note that within each 
layer Xj = x/Dj = kj/kqj where j = 1 or 2) 


( Xp k D j ) = ~ Ei( *DjXj) 


1 p i 

+ - k dj kdj)E\ ( KDj \Xj—Xj I )dxj (7 =1, 2) (5) 

2 J 0 

%{K Dj )=\-2K DJ \ <S>j{Xp K Dj )E 2 (K D jXj)dXj 0 = 1 , 2 ) ( 6 ) 

J o 

It is noted that Eqs. (5) and (6) are the same for each layer, 


Nomenclature 


a 

a , b , c, d 
D 

E\y Ely £ 3 


F(n) 


n 


Q 

q°r 


Tsu 


T 

T s 2 


Xu X 2 


absorption coefficient of 
layer, m “ 1 

the four internal sides of 
the interfaces in the two- 
layer laminate 
thickness of plane layer, 
m 

exponential integral 
functions, 

E m (x) 

= i n m ~ 2 exp(-x/ 'n)dn 
J o 

function of refractive in- 
dex defined in Eq. (15a) 
index of refraction 
heat flux, W/m 2 
externally incident radia- 
tion flux, W/m 2 
absolute temperature, K 
temperatures of sur- 
rounding radiating re- 
gions, K 

coordinates normal to 
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boundaries of individual 
plane layers, m; X\ = 
X\/D\y X 2 = x 2 /D 2 
X = coordinate through two- 
layer laminate, X = X\ 
in first layer, X = 1 + 
X 2 in second layer 
k = optical depth = (a + 
o s )x; k d - optical thick- 
ness = (a + o s )D 
K e = extinction coefficient in 
the complex index of re- 
fraction; n = n - iK e 
o = Stefan-Boltzmann con- 
stant, W/m 2 «K 4 
o s = scattering coefficient of 
layer, m' 1 

p' = reflectivity of interface 
for internally incident 
radiation 

t° - transmissivity of surface 
for externally incident 
radiation 


<t> = dimensionless tempera- 
ture distribution, Eq. (3) 
^ = dimensionless radiative 
heat flux, Eq. (4) 


Subscripts 

iy o = incoming and outgoing 
radiation 

j = index 1 or 2 for each of 
the two sublayers in 
composite 

r = radiative quantity 
1, 2 = layers adjacent to the 
hotter and cooler sur- 
roundings; hotter and 
cooler surroundings 


Superscripts 

/ = at inside surface of an 
interface 

o - at outside surface of an 
interface 

= dummy variable of inte- 
gration 
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and that 4>CA", k d ) and ¥(k d ) are not functions of n . Hence to 
obtain and ^ for all n > 1 it is necessary to solve Eq. (5) 
only once for each k d and use each result to determine ¥ from 
Eq. (6). Tabulated results for some $(X, k d ) and ^(k d ) are 
available from references such as Heaslet and Warming (1965) 
and Siegel and Spuckler (1992). 

Although solutions for $ and ^ are available, the present 
analysis has not yet provided T(ki), T(k 2 ), and q r since Eqs. 
(3) and (4) contain the outgoing boundary fluxes q 0t0 , q Q>b , 
q o c , and q 0td that are unknown. These fluxes must be expressed 
in terms of known quantities; this is accomplished by looking 
at the boundary and interface conditions in detail. 

At the two outer boundaries the internal fluxes are related 
to the transmission of external flux and the reflection of in- 
ternal flux by 

Qo,a~ Qr\Ta + Q\,aPa\ Qo,d = Qr2 7 d^~ Qi y dPd (7tf» b) 

At the inside surfaces of the outer boundaries there are the 
following relations between the radiative flux and the outgoing 
and incoming fluxes: 

Qr = Qo,a - Qi,a\ Qr = “ Qo.d + Qi,d (8a, b ) 


Equations (la) and (8c) are combined to eliminate q iM and 
similarly Eqs. (7b) and (8b) to eliminate q itd . This yields 


Qo,a . / (Qr\ T a QrPa) 

1 “Pa 

(9a) 

Qo,d~ . i (< lr2? d + QrPd ) 

1 ~ Pd 

( 9b) 


By using similar relations at interfaces b and c on either side 
of the internal interface, the following relations are found 
between the outgoing fluxes and the radiant flux q r being trans- 
ferred: 


Qo,b~ . i 1 Qo,c^b~^~ QriPb T’ft)] 
1 ~Pb 

(10a) 

Qo,c~ , i lQo,b r c Qr(Pc r c)l 

1 “ Pc 

(106) 


The radiant flux q r can now be obtained. The q Qta from Eq. 
(9a) and q 0tb from Eq. (10a) are substituted into Eq. (4a). 
The q 0tC is then eliminated by using Eq. (4b), and the q 0td is 
eliminated by using Eq. (9b) . The resulting expression is solved 
for q r to yield 


n o _ , IP- 
i Qr\ . 


?d 


1 -Pa 1 -Pbl-Pd 


7 Qrl 


*1 1 -p‘ Q 1 -Pi, l-p l b \\-p' d * 2 ) 


( 11 ) 


As a result of total internal' reflections at interfaces when 
radiation passes into a medium with a lower refractive index, 
there are the following relations at the four surfaces of the 
interfaces: 


T°a = V-p‘a)n 2 u 




2 


; r2=(l -p' d )n\ 


( \2a-d ) 


Equation (12) is used to eliminate the r’s from Eq. (11) and 
the final result for q r becomes: 



(13a) 

By starting the previous algebra with Eq. (4b) rather than Eq. 
(4a) an alternative expression for q r is found as: 


Qr n\ 

1 , I , ~k , M7 p ° 1 1 

♦2 1-Pc 1 -p‘ d \nj \l-pi *,/ 


(136) 


The q r is the same from expressions (13a) and (13b) ; evaluating 
both relations is a check on the solution. 

The temperature distributions in the two layers are now 
obtained. Starting with Eq. (3b), T 2 (k 2 ) can be obtained from 
n\oT\(K t ) = q oM + (q 0-c - q 0 , d )<t> 2 (* 2 )- The q 0 , c - q oM is 
eliminated by using Eq. (4 b) to give n\aT\{K 2 ) = q od + (q r / 
^ 2 )^ 2 (^ 2 ) - The q 0id is now eliminated by using Eq. (9b), and 
T° d is eliminated by using Eq. (12 d) . The final result for T 2 (k 2 ) 
is placed into the dimensionless form 


o 71(k 2 ) ~ q°2 

Q°r\ - q°rl 


Or Pd ^2(«2> 

q°r\ ~q °2 - Pd 


2 
n 2 


(14 a) 


In a similar fashion the dimensionless T](k\) distribution is 
obtained. Starting with Eq. (3a) and eliminating q 0rb using Eq. 
(4a) gives n\oT\(.K<) = q 0 , a + (<7 r /^i)[4>,(Ki) - 1). The q 0 , a is 
eliminated using Eq. (9a), and T° a is eliminated using Eq. (12a). 
Then -n\q° r2 is added to both sides of the equation. The final 
dimensionless result for T\ (*0 is. 


oT\(K\) q r 2 _ Qr Pa 

Qn-Qri " Q^-Qrill-Pa 


^l(Kl)-l 

*1 


rr 2 


(14 b) 


To use these relations, values of p l are needed for various 
refractive indices. The externally incident radiation is diffuse. 
Although the interfaces are not optically smooth, it is assumed 
that each bit of roughness acts as a smooth facet so that the 
reflectivity can be obtained from the Fresnel interface relations 
for a nonabsorbing dielectric medium. By integrating the re- 
flected energy over all incident directions the relation for p‘(n) 
is (Siegel and Howell, 1981), 




(3* +!)(*-!) 
6(h+1) 2 


+ 


n\n 2 - l) 2 
(/z 2 + l ) 3 



2n\n 2 + 2n-l) 8 kV+1) 

_ _('i 2 + l)(« 4 ‘l) +{n 2 -+l)(n*-l) 2ln{n) - 


(15 a) 


This is for radiation passing into a material of higher refractive 
index where n = n h /n s (n h and n s are the “higher” and 
“smaller” n values). When going in the reverse direction from 
a higher to a smaller n value, the p 1 is given by (Richmond, 
1963) 


P l (n) = 1 — 2 [l-F(n)) ( n = n h /n s ) (15b) 

n 


The relations in Eq. (15) are for perfect dielectrics, that is, 
for materials that do not attenuate radiation internally. As 
discussed by Cox (1965), in the spectral regions where ceramic 
materials are semitransparent to radiation the extinction coef- 
ficient, k € , in the complex index of refraction is usually not 
large enough to significantly affect the surface reflectivity. The 
absorption coefficient, a, in a material is related to K e by a - 
47TK e /Xo. Since wavelengths Xo for thermal radiation are in the 
micrometer range, only a small value of K e is required for this 
relation to yield a large value for a. If the extinction coefficient 
K e is large enough to influence the interface reflectivity rela- 
tions, the absorption coefficient will be so large that the ra- 
diating layer is essentially opaque unless its thickness is much 
smaller than the ceramic layers considered here. 

The required relations have now been provided to evaluate 
Eqs. (13) and (14) for the heat flux and temperature distribution 
in a laminated two-layer region. The equations contain the 
functions $(X, k d ) and V(k d ) that are available in the literature 
as tabulated results for a single layer with n - 1 . For a lam- 
inated layer, these functions are inserted corresponding to the 
optical thickness of each sublayer of the composite. It is not 
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necessary to solve any integral equations to evaluate the present 
solution. 

A convenient alternative to using the exact numerical so- 
lutions for the functions $ and in Eqs. (13) and (14) is to 
use $ and available from the diffusion approximation. The 
accuracy of this approximation is discussed in the next section. 
The diffusion relations for a plane layer are (Siegel and Howell, 
1981) 


*(«/>) = 


*(x, K D )=n*D) 


1 

\ Kl > +X 

1k d (1-X) + 1 - 


(16 a) 


(16 b) 


When used with Eqs. (13) and (14), Eqs. (16) provide rapid 
predictions of heat flux and temperature distributions for any 
optical thicknesses and refractive indices of the sublayers in 
the laminate. 


Table 1 Comparison of results from exact (E) and diffusion (D) solutions 

(a) Dimensionless temperature distribution, <t> 


X 



Optical 

Thickness 



k d =0 

.1 

*0 

= 1 

K 0 = 

10 

K 0 = 

100 


E 

D 

E 

D 

E 

D 

E 

D 

0 

.5710 

.5349 

.7582 

.7143 

.9495 

.9412 

.9948 

.9934 

0.1 

.5541 

.5279 

.6946 

.6714 

.8511 

.8529 

.8974 

.8947 

0.2 

.5397 

.5209 

.6429 

.6286 

.7627 

.7647 

.7994 

.7961 

0.3 

.5262 

.5140 

.5942 

.5857 

.6750 

.6765 

.6998 

.6974 

0.4 

.5130 

.5070 

.5468 

.5429 

.5874 

.5882 

.6000 

.5987 

0.5 

.5000 

.5000 

.5000 

.5000 

.5000 

.5000 

.5000 

.5000 

0.6 

.4870 

.4930 

.4532 

.4571 

.4125 

.4118 

.4000 

.4013 

0.7 

.4738 

.4861 

.4058 

.4143 

.3249 

.3235 

.3002 

.3026 

0.8 

.4603 

.4791 

.3571 

.3714 

.2372 

.2353 

.2007 

.2040 

0.9 

.4459 

.4721 

.3054 

.3286 

. 1488 

.1471 

.1026 

.1053 

1.0 

.4290 

.4651 

.2419 

.2857 

.0505 

.0588 

.0052 

.0066 


(b) Dimensionless heat flux, ¥ 


Results and Discussion 

For a plane layer in radiative equilibrium with absorption 
and isotropic scattering and for n = 1, the $(X, k d ) and ¥(k d ) 
are known functions or can be calculated from Eqs. (5) and 
(6). These functions are the building blocks required to evaluate 
the present solution for a two-layer region including refractive 
index effects and interface reflections. The exact solutions for 
$(X, k d ) and ¥(*p) obtained by solving Eqs. (5) and (6) were 
examined in detail by Heaslet and Warming (1965). The nu- 
merical results used here for $(X, k d ) and ¥(*£>) (Siegel and 
Spuckler, 1992) are in excellent agreement with previous work 
and were carried out to obtain values for a large range of 
optical thicknesses, 0.1 < k d < 100. 

An alternative set of basic functions for $(X, k d ) and V(k d ), 
which was found to be accurate for the present conditions and 
is very simple to implement, is to use the diffusion results in 
Eqs. (16a, b). A comparison of exact and diffusion functions 
is given in Table 1 for various optical thicknesses. The diffusion 
results are reasonably accurate, but there are some differences 
when the optical thickness is small. Results for two-layer lam- 
inated regions were calculated from the present analysis using 
both exact and diffusion functions for the required ${X, k d ) 
and ¥(k d ). There was little difference in the results for the 
laminate, so for convenience the diffusion functions were used 
for evaluating the illustrative results presented here. For sim- 
plicity and sufficient accuracy the diffusion functions are rec- 
ommended for temperature distribution and heat flux 
calculations using the analytical Eqs. (13) and (14). 

The diffusion functions also have an advantage in the theory 
used for developing the present equations for two layers where 
the interfaces are assumed diffuse. Since in the diffusion ap- 
proximation the radiation is everywhere diffuse, there is no 
discontinuity in radiative behavior when energy crosses a dif- 
fuse interface. In the exact solution for a single plane layer 
the diffuse energy leaving an internal boundary is no longer 
diffuse after it travels across the layer. The radiative energy 
passing across from one layer into another thus undergoes a 
change in its angular distribution when crossing a diffuse in- 
terface between the two regions. This can have a small effect 
on the temperatures adjacent to the internal interface. This 
effect is not present when using the diffusion functions for 4> 
and 

For the graphic results shown here, $(X y k d ) and ^( k d ) were 
obtained from Eqs. (16), but exact results such as those in 
Table 1 could have been used. The ^ and were substituted 
into Eq. (13c) to obtain the dimensionless radiative flux through 
the two layers; the heat flux was checked by evaluating the 
alternative Eq. (13d). The dimensionless flux was substituted 
into Eqs. (14 o, b) to obtain dimensionless temperature dis- 


Optical 

Thickness 

** 


E 

D 

0.1 

.9157 

.9302 

1 

.5534 

.5714 

10 

.1168 

.1177 

100 

.0122 

.0132 


Table 2 Interface reflectivities 




P - 1 

Pb 1 

P , 1 

Pd 1 

1 

1 

0 

0 

0 

0 

1 

1.5 

0 

,0918 

.5964 

.5964 

1 

2 

0 

.1606 

.7902 

.7902 

1 

3 

0 

.2762 

.9196 

.9196 

1.5 

3 

.5964 

.1606 

.7902 

.9196 

2 

1 

.7902 

.7902 

.1606 

0 

2 

4 

.7902 

.1606 

.7902 

.9604 

3 

1.5 

.9196 

.7902 

.1606 

.5964 

4 

2 

.9604 

.7902 

.1606 

.7902 


tributions in each layer. The values of the interface reflectivities 
were evaluated from Eqs. (15a) and (15b) depending on 
whether radiation is passing into a region of higher or lower 
refractive index. The reflectivity values at the four surfaces in 
Fig. 1 are given in Table 2 for the combinations of n { and n 2 
values used in Figs. 2-5. In the results that follow, the tem- 
perature distributions are plotted against a continuous X co- 
ordinate where X = X l in layer 1 and X = 1 + X 2 in layer 
2; thus 0 < X < 2. 

Reciprocal Behavior of Dimensionless Tempera- 
tures. Figure 2 illustrates a reciprocal behavior for the tem- 
perature results. Compared to Fig. 2(a ) , both the n values 
and the optical thicknesses of the two layers have been inter- 
changed in Fig. 2(b). It is noted that rotating Fig. 2(a) 180 
deg gives dimensionless temperature distributions that are the 
same as in Fig. 2(b). It was also found that the radiative heat 
flux through the laminate does not depend on the order of 
placement of the two layers relative to the side with larger 
incident radiation. 

Effect of Optical Thickness. Figure 2(a) also illustrates 
the effect of the optical thicknesses of the two layers. The solid 
and dashed profiles are respectively for optical thicknesses of 
the first layer of k d j = 1 and 10. The families of curves cor- 
respond to optical thicknesses in the second layer of k d 2 = 
0.1, 1, 3, 10, 30, and 100. As optical thickness increases, the 
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Fig. 2(a) Layer with lower index of refraction adjacent to higher source Fig. 3(a) Indices of refraction: n, = 1 and n 2 = 2 

of radiant energy 


<D 

3 

"o 

(U 

a 

E 

0) 

£0 

</) 

a; 

c 

o 

’tn 

c 

E 

5 



Fig. 2(b) Demonstration of reciprocity by interchanging n and k d values 


Fig. 2 Dimensionless temperature distributions for indices of refrac- 
tion n, = 1.5 and = 3, two optical thicknesses in layer 1, and six 
optical thicknesses in layer 2 (values reversed in Fig. 2(b)) 



Fig. 3 (b) Indices of refraction: n, = 2 and = 4 


Fig. 3 Effect of magnitude of index of refraction on dimensionless 
temperature distributions in layers as a function of optical thicknesses; 
ratio nj/n, = 2 


resistance to radiative heat flow increases so the temperature 
gradients in the second layer increase. For K m - 0.1 the second 
layer is optically thin and its temperature distribution is es- 
sentially uniform. For k d 2 = 100, however, there is a steep 
gradient in the second layer. When k dx = 10 the large tem- 
perature changes in the first layer cause the temperatures in 
the second layer to be lower than those for k dx = 1. The 
dimensionless radiative heat fluxes are given in Table 3 and 
they illustrate how the heat flux decreases with increasing op- 
tical thickness. 

Effect of Magnitude of Refractive Index. The effect of the 
size of the refractive indices n x and n 2 is illustrated in Figs. 
2(a ) , 3 (a, 6), which form a set of three arrays of temperature 
distributions. The optical thicknesses are as described in the 
previous section. The ratio of n 2 /n x is kept equal to 2 for these 
three sets of profiles, but the refractive indices increase by 
factors of 1:1. 5:2 in the order of Figs. 3(a), 2(a), and 3(b). 
Larger refractive indices increase the amount of internal re- 


flection at the two outer interfaces as a result of total reflection 
when radiation is leaving a material with a higher refractive 
index. This tends to equalize the temperature distributions 
within the layers as the n values increase. The reflections at 
the interface between the layers remain the same because 
n 2 /n x is constant. 

^ For k d i= 1, increasing the refractive index provides de- 
creased temperatures in the first layer. An increased amount 
of the incident radiant energy is reflected away at the first 
interface and does not heat the laminate. Internal reflections 
increase energy absorption within the laminated region, thereby 
increasing the average temperature level in the second layer. 

Effect of Order of Refractive Indices in Heat Flow Direc- 
tion. In Fig. 4 the k d values are the same as in Fig. 3. The 
n values, however, are interchanged so the largest refractive 
index is now adjacent to the environment that provides the 
largest incident radiation. This changes the amount of totally 
reflected energy within the layers as compared with Fig. 3. 
Since total internal reflection is increased in the first layer and 
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Fig. 4(a) Indices of refraction: n, = 2 and /I 2 = 1 



Fig. 4(6) Indices of refraction: n, = 4 and n 2 = 2 


Fig. 4 Effect of inverting the indices of refraction as compared with 
Fig. 3; ratio = 1/2 

decreased in the second, the temperatures become more uni- 
form in the first layer and have steeper gradients in the second. 
The temperature profiles in the first layer are not as close 
together as on Fig. 3. 

Effect of Changing Refractive Index in One Layer. For all 
the distributions in Fig. 5 the refractive index in the first layer 
is retained at n x - 1. In the second layer the values are: n 2 = 
1, 1.5, 2, and 3. The optical thicknesses in the two layers are 
equal. Three families of curves are shown for k d 1 = kq 2 = 
0.2, 2, 20. Hence the solid curves for each layer are optically 
thin, the dashed curves are for an intermediate optical thick- 
ness, and the short dashed curves are for a rather thick region. 

Consider the solid profiles for k D \ = k&i - 0.2. When n = 
1 there are no interface reflections and a continuous temper- 
ature profile is obtained through the laminated region. For n 2 
> 1 the profiles become discontinuous at the internal interface, 
and as n 2 becomes larger the discontinuity increases. For k D \ 
- k d 2 = 0.2, an increased n 2 raises the temperatures in the 
first layer and makes the temperatures more uniform in the 
second layer. For the short dashed curves with k D i - k d 2 = 
20, there is a different effect. As n 2 increases, the temperatures 



Dimensionless coordinate, X 

Fig. 5 Effect of index of refraction ratio on dimensionless tem- 
perature distributions for optical thickness varying from thin to thick 


Table 3 Dimensionless radiative heat flux, ?,/(<£ - <£) 




= 

1 



<02 

n, = 1 
_ n; = 2 

n, = 1.5 
_ n 2 ._= 3 

n, = 2 
n ? = 4 

n i = 2 

n ? = 1 

nj =■ 4 
n * = 2 

0.1 

0.3447 

0.3574 

0.3393 

0.4175 

0.3545 

1 

0.3257 

0.3480 

0.3345 

0.3257 

0.3345 

3 

0.2903 

0.3290 

0.3243 

0.2108 

0.2972 

10 

0.2102 

0.2760 

0.2931 

0.1018 

0.2138 

30 

0.1175 

0.1890 

0.2299 

0.0403 

0.1187 

100 

0.0462 

0.0899 

0.1311 

0.0129 

0.0464 




K oi_ 

= 10 



*32 

ni = 1 
n, = 2 

n : = 1.5 
n, = 3 

n, = 2 
n, = 4 

n i = 2 

n, = 1 

nj = 4 
n 7 = 2 

0.1 

0.1036 

0.1725 

0.2158 

0.2449 

0.3084 

1 

0.1010 

0.1703 

0.2138 

0.2102 

0.2931 

3 

0.0981 

0.1656 

0.2096 

0.1598 

0.2641 

10 

0.0869 

0.1510 

0.1961 

0.0869 

0.1961 

30 

0.0655 

0.1206 

0.1657 

0.0377 

0.1130 

100 

0.0352 

0.0708 

0.1073 

0.0127 

0.0455 


_ », = 1 

_ n, 

A 

11 

a* 

II 

O 

Km = *„ = 2 

= K„ = 20 

1 

0.7692 

0.2500 

0.0323 

1.5 

0.5065 

0.2548 

0.0427 

2 

0.4310 

0.2495 

0.0479 

3 

0.3548 

0.2316 

0.0518 


in the first layer are reduced, while the temperatures in the 
second layer become more uniform as before. For an inter- 
mediate optical thickness, k di - k D2 = 2, there is a small 
decrease in temperatures in the first layer as n 2 is increased, 
and then with a further increase in n 2 the temperatures increase. 
The radiative heat fluxes in Table 3 show a similar reversal in 
trend when going from an optically thin to an optically thick 
layer. 

Conclusions 

A convenient solution was obtained for temperature distri- 
butions and radiative heat fluxes in a laminated two-layer 
region with radiative energy incident on each outer boundary. 
The interfaces are assumed diffuse with the intent of modeling 
interface reflections of a ceramic material. The limiting case 
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considered here is for a two-layer laminate in radiative equi- 
librium (no heat conduction) that emits, absorbs, and iso- 
tropically scatters radiation. The solution yields results for two 
laminated layers by using known dimensionless temperature 
and heat flux functions for a single semitransparent radiating 
layer with n - 1. These functions are combined through the 
coupling at the interfaces to yield the solution for a two-layer 
laminated region with n > 1 . The coupling relations include 
transmission through the interfaces and all internal reflections 
within the layers. The result is a set of simple algebraic expres- 
sions for the temperature distribution in each of the layers, 
and for the radiative heat flux through the laminated region. 
Results were evaluated using the basic functions for a single 
layer obtained by both exact and diffusion methods. Unless 
the optical thickness is much smaller than one, the diffusion 
results are within engineering accuracy, and they provide an 
especially simple way to compute results for two layers. 

Illustrative results are given for a variety of refractive indices 
and optical thicknesses of the layers. A reciprocity relation 
was found for the temperature distributions. If the order of 
the refractive indices of the layers is reversed as well as their 
optical thicknesses, the temperature distributions become in- 
verted mirror images. It was found that increasing the refrac- 
tive index in each layer, while keeping their ratio constant, 
generally makes the temperature profiles more uniform as the 
result of increased internal reflections. When the refractive 
index of the second layer is increased, and the first is kept 
unchanged, the temperatures in the layers can either increase 
or decrease depending on the layer optical thicknesses. 
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